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Point defect complexes in crystalline silicon composed of hydrogen, nitrogen, and oxygen atoms are studied 
within density-functional theory (DFT). Ab initio Random Structure Searching (AIRSS) is used to find low- 
energy defect structures. We find new lowest-energy structures for several defects: the triple-oxygen defect, 
{30i}, triple oxygen with a nitrogen atom, {Ni,30i}, triple nitrogen with an oxygen atom, {3Ni,Oi}, double 
hydrogen and an oxygen atom, {2Hi,Oi}, double hydrogen and oxygen atoms, {2Hi,20i} and four hydro- 
gen/nitrogen/oxygen complexes, {Hi, Ni, Oj}, {2Hi, Oj}, {Hi, 2Nj, Oj}, and {H i5 Ni, 20i}. We find that some 
defects form analogous structures when an oxygen atom is replaced by a NH group, for example, {Hi,Ni,20i} 
and {30i}, and {Hi,Ni} and {Oi}. We compare defect formation energies obtained using different oxygen 
chemical potentials and investigate the relative abundances of the defects. 

PACS numbers: 61.05.-a, 61.72.jj 



I. INTRODUCTION 

Hydrogen (H) is a common impurity in silicon (Si) which 
is particularly important as a passivator of surfaces and bulk 
defects. 1 The role of hydrogen in semiconductors is high- 
lighted in a review by Estreicher. 2 Adding nitrogen (N) im- 
purities to silicon affects the formation of voids and may in- 
crease the strength of the silicon crystal by immobilising dis- 
locations, which reduces warping during wafer processing. 3 
The majority of silicon used in device technologies is manu- 
factured in a quartz crucible by the Czochralski (Cz) process, 
during which oxygen (O) from the quartz readily enters the 
melt, see for example the review by Newman. 4 

A wide variety of experimental probes are used to study 
defects in silicon, but it is often difficult to determine their de- 
tailed structures from measurements alone. Theoretical stud- 
ies using ab initio techniques are helpful in this regard, as 
they are used to calculate both the defect formation energies 
and some of their experimental signatures. For example, local 
vibrational modes of defects are accessible to infra-red (IR) 
absorption experiments and may also be calculated within ab 
initio methods. 

Throughout this paper we denote a defect complex by list- 
ing its constituent atoms between braces, {Ai,# s , . . .}, where 
a subscript denotes whether an atom is substitutional (s) or in- 
terstitial (i). For example, a defect containing two interstitial 
hydrogen atoms and a substitutional nitrogen atom is denoted 
by {2Hi,N s }. Despite hydrogen and nitrogen being common 
impurities in silicon, we are aware of only one ab initio the- 
oretical study of their interaction. 5 It was concluded that both 
{N s } and {Ni} defects are able to trap hydrogen atoms, al- 
though with smaller binding energies than that of a hydrogen 
atom and a vacancy. The Fourier transform infra-red (FTIR) 
vibrational line at 2967 cm -1 found in silicon prepared in a 
hydrogen atmosphere is assigned to the N-H stretch mode of 
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the {Hi,Ni} defect. 6 Interstitial oxygen may be present in Cz- 
Si in concentrations as high as 10 18 cm -3 . Oxygen impurity 
atoms may be used for gettering metallic impurities, which 
increases the overall crystal quality, but they also form elec- 
trically active thermal donors (TD). 7 TDs can affect the local 
resistivity within the silicon wafer. 8 In addition, oxygen also 
interacts with nitrogen in bulk silicon to form TDs, and small 
concentrations of hydrogen can greatly enhance the forma- 
tion of TDs. 9 It has been shown that the presence of hydrogen 
molecules in silicon crystals enhances the diffusion of oxy- 
gen. It has also been suggested that hydrogen passivates the 
electrical activity of N/O complexes in Cz-Si. 10 

It is clear from these examples that the interactions between 
hydrogen, nitrogen, and oxygen impurities in silicon lead to 
important and complicated behavior. Here we present the re- 
sults of a computational search for low-energy defect com- 
plexes in silicon containing hydrogen, nitrogen, and oxygen 
atoms. We confirm the stability of many of the previously- 
known defects and also report the structures and formation 
energies of some new low-energy defects. To the best of our 
knowledge, this is the first time that all of these defects have 
been compared within a single study. 

The rest of this paper is set out as follows, in Sec. (TTJ we 
discuss the computational methods used in the study. In Sec. 
Hill we describe the calculation of the defect formation ener- 
gies and explain our choice of chemical potentials. Results 
for H/N defects are presented in Sec. [iVl N/O defects in Sec. 
El H/O defects in Sec.ED and H/N/O defects in SecED The 
relative abundances of the defects are studied in Sec. I Villi and 
a discussion of our results is presented in Sec.HXl 



II. COMPUTATIONAL APPROACH 

A. Random Structure Searching 

"Ab initio random structure searching" (AIRSS) has al- 
ready proved to be a powerful tool for finding crystalline 
structures of solids under high pressures ^MWW 6 The 
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basic algorithm is very simple: we take a population of ran- 
dom structures and relax them to local minima in the energy. 
Pickard and Needs 11 showed that "random structure search- 
ing" can be applied to finding structures of point defects, and 
the approach is discussed in more detail in Ref . |_1V|. 

Creating the initial simulation cell is a three-stage process. 
First, we choose the size of the disruption to the perfect silicon 
lattice by defining the radius of a sphere in which the impurity 
atoms are to be placed randomly. We choose radii between 3 
and 7 A. Secondly, we chose the number of silicon atoms to 
have their positions ramdomized. In this study we randomize 
the position of 1 to 3 silicon atoms depending on the num- 
ber of hydrogen, nitrogen and oxygen impurity atoms present. 
The randomization sphere is centered on the centroid of the 
removed silicon atoms. Finally, the required number of hy- 
drogen, nitrogen and oxygen impurity atoms and the chosen 
1 , 2 or 3 silicon atoms are placed at random positions within 
the sphere. The configurations are then relaxed using the cal- 
culated density-functional-theory (DFT) forces. 



B. Density-functional-theory calculations 

Our calculations are performed using the General- 
ized Gradient Approximation (GGA) density functional 
of Perdew, Burke, and Ernzerhof (PBE). 18 We use the 
plane- wave basis-set code CASTEP, 19 with its built-in ul- 
trasoft pseudopotentials 20 which include non-linear core 
corrections. 21 All of the results presented here are for 
non- spin-polarized systems. The multi-Baldereschi k-point 
scheme described in Ref.[l7|is used for sampling the Brillouin 
zone. We perform the same computational convergence tests 
as described in Ref. [l7l finding excellent convergence with 
a 2x2x2 k-point sampling, which we use for all the results 
reported here. 

The searches we report in this paper are carried out in 
a body-centered-cubic simulation cell which would hold 32 
atoms of crystalline silicon. Although this is a small cell, 
tests indicate that it is adequate for obtaining good structures 
for most of the defects studied and reasonably accurate for- 
mation energies, while keeping the computational cost low 
enough to permit extensive searching. In a previous study, 17 
we reported results for hydrogen/silicon complexes in silicon. 
The searches were carried out in the same 32-atom cell used 
here, but each of the defect structures was re-relaxed in a 128- 
atom cell to obtain more accurate structures and energies. The 
defect formation energies for the 32-atom cells were not re- 
ported, but they are in fact in excellent agreement with the 
128-atom values. Consider, for example, the calculations for 
an interstitial silicon atom (denoted by /) and hydrogen im- 
purity atoms in crystalline silicon reported in Ref. 17. The 
changes in the formation energies on increasing the cell size 
from 32 to 128 atoms for the defects {/,H}, {/,H 2 }, {/,H 3 } 
and {/,H 4 } are only 0.00, -0.07, -0.01 and -0.02 eV, re- 
spectively. We have also compared our results for N/O defects 
in silicon with data from the literature. Fujita et al^ used 
DFT and 216-atom cells to calculate the binding energies for 
the following defect reactions in silicon: {2Ni,Oi}+{Oi} — > 



{2Ni,20i}, {2Ni}+{Oi} -> {2Ni,Oi}, and {NiA} + {Oi} 
-> {Ni,20i}, obtaining 0.91, 0.98, and 1.12 eV, respectively, 
in good agreement with our 32-atom cell results of 0.80, 0.91 
and 0.97 eV, respectively. Finally, Coutinho et al. 23 used DFT 
and 64-atom cells to calculate the formation energy of the 
{Oi} defect, obtaining values between 1.989 and 1.820eV, 
depending on the basis set used, which is in good agreement 
with our value of 1.90 eV. We should not, of course, expect a 
32-atom cell to be sufficient for all point-defect calculations in 
silicon, and especially not for charged defect states, but there 
is a strong case for believing it to be adequate for the purposes 
of the current study. 



III. CALCULATING THE FORMATION ENERGIES 

In order to define defect formation energies it is necessary 
to specify the chemical potentials of the atomic species. The 
chemical potential for silicon is taken to be the bulk silicon 
value, so that the energy cost to add a silicon atom to a de- 
fect is the energy of a silicon atom in perfect crystalline sili- 
con. The chemical potential for hydrogen is derived from the 
lowest energy hydrogen defect in silicon, which is an inter- 
stitial molecule at the tetrahedral site, {H2 i}, with the bond 
pointing along a (100) direction. This is the same chemi- 
cal potential as used in our study of hydrogen impurities in 
silicon. 17 We use the most energetically favorable nitrogen 
defect, {2Ni}, 5 ^ 25 ^ 27 to set the chemical potential of ni- 
trogen. 

We perform around 130 AIRSS searches with one, two, 
three, or four interstitial oxygen atoms per cell. The lowest en- 
ergy {Oi} defect found is the buckled-bond-centered config- 
uration which has been studied extensively. 23,28,29,30,31,32 The 
{20i} defect has two buckled-bond-centered oxygen atoms 
bonded to a silicon atom. This defect was reported in Ref. HI 
and is known as the staggered {20i} defect. 

Our searches also find the previously-reported staggered 
chain {30i} defect, 30 ^ 3 ^ 4 ^ 5 ^^ 37 We find a new structure, see 
Fig. [U with a formation energy 0.15 eV lower than the stag- 
gered chain. To the best of our knowledge this defect has not 
been mentioned before. We also find the {40i} defect which 
was previously reported as the 04(1,1) defect in Ref. [38l. 

Table U shows the formation energies per oxygen atom of 
the four oxygen defects and quartz, using the chemical poten- 
tial of the {Oi} defect as the reference for supplying oxygen 
atoms. We note that the formation energy decreases with each 
additional oxygen, showing that oxygen defects tend to aggre- 
gate, the lowest possible energy being achieved by forming 
crystalline Si02 (quartz). 

We consider three different chemical potentials for oxygen. 
These correspond to choosing the source of oxygen atoms as 
/i{o-}, M{40/}» an d Mquartz- The chemical potential jU{ n o/} f° r 
oxygen with the source of oxygen atoms being the {nO{} de- 
fect is calculated as the energy of the lowest-energy struc- 
ture of the 32-atom silicon cell containing n oxygen atoms 
and the energy of the 32-atom bulk silicon cell, divided by n. 
The {Oi} defect gives the highest oxygen chemical potential, 
while quartz gives the lowest chemical potential. We include 
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the chemical potential, /i qua rtz, calculated from quartz, since it 
is the opposite extreme to /i{o, } m terms of oxygen saturation. 



Defect Ef per O (eV) 
O 0.00 
20 -0.25 
30 -0.42 
40 -0.46 
Quartz -1.90 



TABLE I: Formation energies per oxygen atom for different oxygen 
complexes, relative to the oxygen bond-centered defect (Ot, c )- These 
values can also be thought of as the oxygen chemical potentials rela- 
tive to a source of oxygen atoms consisting of single-oxygen defects. 




FIG. 1: (Color online) The {30i} defect is composed of three 
buckled-bond-centered oxygen atoms, each bonded to a single, four- 
fold coordinated silicon atom. Silicon atoms are shown in yellow and 
oxygen atoms in red. 
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TABLE II: Formation energies defined by Eq. in eV for vari- 
ous H/N/O complexes with three different choices of the oxygen 
chemical potential. The final column indicates whether the system 
could form a structure with fully- saturated covalent bonding (y0> or 
whether it cannot (x), see Sec. IIXI The defect structures shown in 
bold have not, to the best of our knowledge, been reported in the 
literature before. 

The formation energy of a defect, Ef , is defined as 

E f = E D -Y^niHi, (1) 

i 

where U[ is the number of each atomic type i in the defect 
cell of energy and the chemical potentials, are defined 
above. 




FIG. 2: (Color online) The first metastable H/N defect, {2Hi,N s }, 
contains a three-fold coordinated nitrogen atom close to a lattice site. 
One of the hydrogen atoms is bonded directly to the nitrogen atom 
and points toward a neighboring three-fold coordinated silicon atom. 
The remaining hydrogen atom terminates a dangling bond on another 
neighboring silicon atom. The silicon atoms are shown in yellow, the 
nitrogen in blue and the hydrogen atoms in white. 




FIG. 3: (Color online) The {Ni, 30i} defect contains a four- atom ring 
consisting of a nitrogen atom, an oxygen atom and two silicon atoms. 
The two other oxygen atoms are in buckled-bond-centered positions. 
Note that the structure contains an over-coordinated oxygen atom 
with three bonds. The silicon atoms are shown in yellow, the nitrogen 
in blue, and the oxygen atoms in red. 
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FIG. 4: (Color online) The lowest-energy {3Ni,Oi} defect consists 
of a four- atom ring of two nitrogen atoms and two silicon atoms 
adjacent to another four- atom ring of an oxygen atom, a nitrogen 
atom and two silicon atoms. Note that the structure contains an over- 
coordinated oxygen atoms with three bonds. The silicon atoms are 
in yellow, the nitrogen atoms in blue and the oxygen atom in red. 




FIG. 5: (Color online) The {2Hi,Oi} defect is composed of a 
buckled-bond-centered oxygen atom adjacent to a distorted {H2}* 
defect. The silicon atoms are shown in yellow, the oxygen atom in 
red and the hydrogen atoms in white. This defect has a negative 
formation energy, showing that the metastable {H2}* defect 39 is sta- 
bilized by the presence of oxygen. 




FIG. 6: (Color online) The {2Hi,20i} defect is composed of two 
adjacent {Hi,Oi} defects. The silicon atoms are shown in yellow, 
the oxygen atom in red and the hydrogen atoms in white. 




FIG. 7: (Color online) The {Hi,Ni,Oi} defect contains a nitrogen 
and silicon atom sharing a lattice site, with the dangling silicon bond 
terminated by the hydrogen atom. A nearest-neighbor silicon atom 
of the nitrogen atom is bonded to a buckled-bond-centered oxygen 
atom. The silicon atoms are shown in yellow, the oxygen in red, the 
nitrogen in blue and the hydrogen in white. 




FIG. 8: (Color online) The {2Hi,Ni,Oi} defect is similar to the 
{Hi,Ni,Oi} defect but with the nitrogen atom acquiring the extra 
hydrogen atom, making it over-coordinated with four bonds. This 
defect has the highest formation energy of all the H/N/O defects stud- 
ied and is therefore unlikely to form. The silicon atoms are shown 
in yellow, the oxygen in red, the nitrogen in blue and the hydrogen 
atoms in white. 



IV. H/N DEFECTS 

We have relaxed around 250 different starting structures for 
the H/N defects using cells containing 32 silicon atoms, 1 ni- 
trogen atom and 1 hydrogen atom, and further searches with 
31 silicon atoms, 1 nitrogen atom and 1 hydrogen atom, and 
with 31 silicon atoms, 1 nitrogen atom and 2 hydrogen atoms. 
These are the only calculations reported in this paper where 
we change the number of silicon atoms in the cell. We find 
defects which we may describe as bulk silicon with an intersti- 
tial nitrogen atom and an interstitial hydrogen atom, {Hi,Ni}, 
bulk silicon with a substitutional nitrogen atom and an inter- 
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FIG. 9: (Color online) The {Hi, 2Ni, Oi} defect contains a four-atom 
ring of two nitrogen atoms and two silicon atoms. One of the nitrogen 
atoms is over-coordinated and is bonded to the hydrogen atom. A 
nearest-neighbor silicon atom to this nitrogen atom is also bonded to 
a buckled-bond-centered oxygen atom. The silicon atoms are shown 
in yellow, the oxygen in red, the nitrogen in blue and the hydrogen 
in white. 




FIG. 10: (Color online) The {Hi,Ni,20i} defect is very similar to 
the {30i} defect, but with one of the oxygen atoms replaced by a 
NH group. The silicon atoms are shown in yellow, the oxygen in red, 
the nitrogen in blue and the hydrogen in white. 



stitial hydrogen atom, {Hi,N s }, and bulk silicon with a sub- 
stitutional nitrogen atom and two interstitial hydrogen atoms, 
{2Hi,N s }. 

The most stable H/N structure we find is the {Hi,Ni} de- 
fect with C s symmetry. In this defect the hydrogen atom is 
bonded to the buckled-bond-centered nitrogen atom, which is 
three-fold coordinated, the nitrogen, hydrogen, and two sili- 
con atoms lie in a plane. The N-Si bond lengths are both about 
1 .7 1 A and the N-H bond-length is 1 .03 A, in good agreement 
with McAfee et al., 5 who found bond lengths of 1.73 A and 
1.05 A, respectively. 

The first metastable defect we find, {2Hi,N s }, is shown in 
Fig. [2] and has C^ v symmetry. This is quite a complicated de- 
fect, and it is energetically slightly more favorable than the 



{Hi,N s } defect which is described below. {2Hi,N s } looks 
similar to the {Hi, N s } defect but with the extra hydrogen atom 
bonded to the three-fold coordinated nitrogen atom. 

Finally, we report the {Hi,N s } defect. This is similar to 
the C3 V symmetry {N s } defect, 5 but with the hydrogen atom 
bonded to the three-fold coordinated silicon atom, increasing 
its coordination number to four. The three N-Si bonds are all 
1.85 A long, the same as in {N s }. The N-H distance is 1.95 A 
and the H-Si bond is 1.80 A long, in good agreement with 
McAfee et al 5 

All of the H/N complexes we mention in this section have 
Ef > 0, and hence they are unlikely to form spontaneously. 



V. N/O DEFECTS 

Our searches for the N/O complexes required around 
1800 random starting structures. These searches found all 
the previously-known lowest-energy structures except the 
{2Ni,20i} defect. The {2Ni,20i} defect is quite large and it 
is probable that we have not relaxed a sufficiently large num- 
ber of starting structures to find it. 

For the {Ni, Oi} defect we obtain the previously-known in- 
terstitial ring. 40 This defect has a positive £f(O0 of 0.43 eV, 
and hence is unlikely to form. The structure of our {Ni,20i} 
defect is the same as found in previous studies. 22 i 4Q i 41 This de- 
fect has a negative Ef(0[) of — 0.54eV, and therefore it could 
form in silicon. 

We are not aware of any previous reports of the {Ni,30i} 
defect, which is shown in Fig. [3j The {Ni,30i} defect has 
a negative Ef(0{) of — 0.92eV, showing that a single bond- 
centered oxygen atom will bind to the {Ni, 20i} defect. 

Our searches for the {2Ni,20i} defect did not yield the 
structure reported previously, 42 which has different oxygen 
positions. We generated the structure reported by Fujita et 
al 42 "by hand" and found it to be 0. 19 eV lower in energy than 
our best structure. The defect of Fujita et al— has a negative 
formation energy of Ef(0{) = — 1 .71 eV, and the formation en- 
ergy remains negative even when the reference structure for 
the oxygen chemical potential is taken to be that of the {40i} 
defect in silicon. It is disappointing that our searches has not 
given the previously-known lowest-energy {2Ni,20i} defect. 
However, we have included the lowest-energy known struc- 
ture in our analysis of the relative populations of the various 
defects presented in Sec. I Villi 

We find the previously-known lowest-energy {2Ni,Oi} 
defect j 4( W 42 i 43 which has a formation energy of £f(Oi) = 
—0.91 eV. This defect binds oxygen very strongly and £f(Oi) 
remains negative even with the oxygen chemical potential 
from the {40i} defect. 

We are not aware of any previous reports of the {3Ni,Oi} 
defect in the literature, which is shown in Fig. [4j The 
{3Ni, Oi} defect is quite large and it may not be well described 
within a 32-atom cell. 
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VI. H/O DEFECTS 

We have relaxed around 350 starting structures for the H/O 
defects. The lowest energy {Hi,Oi} defect that we find is 
composed of a buckled-bond-centered oxygen atom adjacent 
to a buckled-bond-centered hydrogen atom. This structure 
was also found by Estreicher i 44 i 45 This defect has, however, 
been the subject of some controversy, as Jones et al 46A1 pro- 
posed that the hydrogen atom is attached to a silicon atom at 
an anti-bonding site. Our search also finds this defect config- 
uration, but we calculate it to be 0.32 eV higher in energy than 
the ground state structure. 

Our lowest-energy {Hi, Oi} defect has a positive formation 
energy, and hence it is unlikely to form. The {2Hi, 20i} defect 
is composed of two such defects in close proximity to one an- 
other, as shown in Fig. [6] The defect has a negative formation 
energy and, to our knowledge, it has not been presented in the 
literature previously. 

The lowest-energy {2Hi,Oi} defect that we find (Fig. 
adopts a {H2}*+Ob c configuration. Measurements have 
shown an IR absorption line at 1075.1 cm -1 which has been 
assigned to a Oi-H2 complex. 48 However, the {H2}*+Ob c con- 
figuration of the {2Hi, Oi} defect does not contain a hydrogen 
molecule, see Fig. [5] We also find a {2Hi, Oi} defect structure 
containing a hydrogen molecule, but it is metastable with an 
energy 0.26 eV above our ground-state structure. To the best 
of out knowledge the {H2}*+Ob c defect structure has not been 
reported previously. 

Our most stable {Hi,20i} defect is a Ob c +Ob c +H configu- 
ration. This defect has also been studied by Rashkeev et al 
(see Fig. 3(a) of Ref. l36|) . who report it to be a thermal double 
donor (TDD). 

The {Hi,20i}, {2Hi,Oi} and {2Hi,20i} defects all have 
negative formation energies Ef(0{) and are therefore likely to 
form in bulk silicon. 



VII. H/N/O DEFECTS 

We have performed around 500 structural relaxations for 
the H/N/O defects. The four defects {Hi,Ni,Oi} (Fig. 0, 
{2Hi,Ni,Oi} (Fig.©, {Hi,2Ni,Oi} (Fig.0, and {H^N^Oi} 
(Fig. [10]), all have negative formation energies when using 
Ef(0{), and {Hi,Ni,Oi} and {Hi,2Ni,Oi} also have negative 
formation energies when using £f(40i). 

The {Hi, Ni, Oi} defect has a formation energy of —0.68 eV 
compared with 0.43 eV for {Ni, Oi}. This implies that hydro- 
gen readily binds to the {Ni,Oi} defect. The hydrogen atom 
bonds to a silicon atom and breaks the third bond of the over- 
coordinated oxygen atom, giving the structure shown in Fig. 
[71 This results in a defect with fully saturated bonds which is 
therefore quite low in energy. 

The {2Hi,Ni,Oi} defect is very similar to the {Hi,Ni,Oi} 
defect mentioned above. The extra hydrogen atom bonds to 
the nitrogen atom which is then over-coordinated, having four 
bonds, see Fig. [8] The over-coordinated nitrogen atom is en- 
ergetically unfavorable and this defect would not readily form 



from {Hi,Ni,Oi}, however it could form from {Ni,Oi} by 
capturing a hydrogen molecule. 

The {Hi,2Ni,Oi} defect is an interesting case. The nitro- 
gen atom is over-coordinated because it binds to the hydro- 
gen atom, see Fig. [9] Since {2Ni,Oi} has a formation en- 
ergy of — 0.9 leV and {Hi,2Ni,Oi} has a formation energy 
of — 0.49eV it is unlikely that a hydrogen atom will bind to 
{2Ni,Oi}. 

The {Hi,Ni,20i} defect is very interesting. The fully- 
saturated bonding of this defect results in a low formation 
energy of —0.85 eV. This is lower than the formation energy 
of {Ni,20i}, and hydrogen will therefore bind to this defect. 
The {Hi,Ni,20i} defect has a similar geometry to {30i}, but 
with one of the oxygen atoms replaced by a HN group. A 
HN group acting in a similar fashion to an oxygen atom was 
also observed in the {Ni,Hi} defect of McAfee, 5 which looks 
similar to the bond-centered {Oi} defect. 



VIII. RELATIVE ABUNDANCES 

The formation energies of the various defects can be used to 
calculate the relative abundances of the defects at zero temper- 
ature as a function of the ratios of the H/N/O concentrations. 
The relative abundances for a chosen set of concentrations are 
those which minimise the total energy, and the chemical po- 
tentials do not enter the calculation. However, if we allow 
the oxygen atoms to combine with silicon to form quartz, it 
is straightforward to show that the minimum energy is always 
obtained by producing as much quartz as possible, so that no 
point defects containing oxygen atoms remain in the bulk, and 
the hydrogen and nitrogen atoms form {H2 i} and {2Ni} de- 
fects. This is inconsistent with the extensive experimental ev- 
idence for other point defects in silicon involving hydrogen, 
nitrogen, and oxygen atoms. To gain some insight into the 
H/N/O defects which may be present we must therefore limit 
the propensity to form quartz and other low-energy oxygen 
defects in some way. We have chosen to present relative abun- 
dances in which the formation of {20i}, {30i}, {40i}, up 
to quartz are excluded, while allowing the formation of {Oi}. 
This is an arbitrary choice, but it serves to illustrate the type of 
behavior which may arise. Figs. [TT] [321 and [T3l show relative 
abundances of the defects as a function of the H/N/O concen- 
trations. In each figure we keep the concentrations of two of 
the species constant and equal to unity, while the concentra- 
tion of the third impurity varies from zero to three. The main 
features of Figs. [TT] [121 and [13] are that, for equal concentra- 
tions of hydrogen, nitrogen, and oxygen, the dominant defects 
are hydrogen molecules and the {2Ni,20i} defect. Varying 
the nitrogen or oxygen concentration leads to the formation 
of a wide variety of defects. Increasing the hydrogen concen- 
tration, however, only generates more hydrogen molecules. 
Other scenarios can be investigated using the data given in 
Table M 
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FIG. 1 1 : (Color online) Relative concentrations of H/N/O defects 
in silicon as a function of the nitrogen concentration. At low nitro- 
gen concentrations (H:N:0 = 1:<0.5:1), the predominant defect is 
{2Hi,20i}. As the nitrogen concentration increases the concentra- 
tion of {2Hi,20i} decreases, with {2Ni,20i} and molecular hydro- 
gen being formed. As the nitrogen concentration increases above 
unity, the concentration of {2Ni,20i} declines and {2Ni,Oi} begins 
to form. At nitrogen concentrations larger than two, the additional 
nitrogen atoms form {2Ni} defects. 
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FIG. 12: (Color online) Relative concentrations of H/N/O defects in 
silicon as a function of the oxygen concentration. The behavior is 
complex in this case, although the general trend is simply that the 
more stable defects containing larger numbers of oxygen atoms are 
favored at higher oxygen concentrations. 



IX. DISCUSSION 

We have presented ab initio calculations of H/N/O com- 
plexes in bulk silicon within DFT. The defect complexes were 
generated in 32-atom supercells using the AIRSS method as 
outlined in Ref . The 32-atom cells are large enough to 
give a reasonable description of most of the defects, but small 
enough to allow extensive searching. To the best of our knowl- 
edge this it the first time that H/N/O complexes in silicon have 
been studied in detail. 

AIRSS searches were performed on oxygen, hydrogen, and 
nitrogen defects to evaluate chemical potentials for the atoms. 



FIG. 13: (Color online) Relative concentrations of H/N/O defects in 
silicon as a function of the hydrogen concentration. Increasing the 
hydrogen concentration simply generates more interstitial hydrogen 
molecules. 



Whilst carrying out the searches on oxygen defects, we dis- 
covered a new {30i} defect, shown in Fig. [U which is more 
stable than the staggered defect which has been reported in the 
literature. 

Our searches reproduced all of the previously-known 
lowest-energy H/N complexes. We have also described a 
new metastable {2Hi, N s } defect which is more favorable than 
some of the other H/N complexes presented in the literature. 

The searches for N/O complexes were slightly less success- 
ful. We were unable to find the previously-known {2Ni,20i} 
defect. In this case, we performed DFT calculations on the 
previously-known best structure for comparison. We discov- 
ered a new {3Ni,Oi} defect (Fig. [4]) and a new {Ni,30i} de- 
fect (Fig. [3]) which is more stable than the {Ni,20i} defect, 
implying that an oxygen atom could bind to {Ni, 20i}. 

Our searches reproduced all of the previously-known 
lowest-energy H/O complexes. We find the {Hi,Oi} defect 
to have a positive formation energy, and therefore it is un- 
likely to form, but two adjacent {Hi,Oi} defects can form a 
{2Hi,20i} defect (Fig.O which has a formation energy of 
£f(Oi)= -0.76 eV. We also find a new {2Hi, Oi} defect which 
is based on the previously-reported metastable HJ defect. 39 
The {2Hi, Oi} defect has £ f (Oi)= -0.23 eV, showing that the 
HJ defect can bind an oxygen atom. 

To the best of our knowledge, no previous discussion of 
H/N/O complexes in silicon has appeared in the literature. 
The {Hi,Ni,Oi} and {Hi,Ni,20i} defects are more energet- 
ically favorable than their non-hydrogen-containing counter- 
parts {Ni,Oi}, {Ni,20i}, showing that many of these defects 
could exist in hydrogenated forms. 

Another point of interest is that H/N complexes behave in 
an analogous fashion to oxygen defects. For example, the 
{Hi,Ni} defect adopts a buckled-bond-centered arrangement 
with the hydrogen atom terminating the dangling bond on the 
nitrogen atom. The nitrogen atom takes a similar position to 
the oxygen atom in the Ot> c defect, however without the hydro- 
gen atom, since oxygen forms only two covalent bonds. The 
same relationship holds between the {30i} and {Hi,Ni,20i} 
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defects, where again the H/N group in the latter behaves like 
an Obc defect in the former. 

The right-hand column of Table HT1 indicates whether or not 
it is possible for the indicated set of atoms to form a fully- 
saturated covalently-bonded structure. In such a structure 
each silicon atom is bonded to four other atoms, each nitro- 
gen to three, each oxygen to two, and each hydrogen to one. 
As each bond is shared between two atoms, a fully- saturated 
covalently-bonded structure is impossible if 

4N S i + 3N N + 2N + /% = odd number, (2) 

where A/x denotes the number of atoms of type X in the cell. 
We note from the data in Table HH that the lowest energy de- 
fects are normally those having completely saturated bonds. 
A similar result for hydrogen defects in silicon was reported 
in Ref . |_1V|- This simple and rather obvious rule is useful for 
choosing combinations of impurity atoms which might form 
stable defects. 

Consider a host crystal containing three types of impurity 
atom. (An example might be H/N/O impurities in silicon as 
studied here, but without the possibility of silicon vacancy for- 
mation which was allowed in our study of H/N defects.) Lim- 
iting the total number of impurity atoms to be <4 gives 34 
possible cell contents for which searches must be performed. 
Increasing the number of possible impurities to five (while 
maintaining the maximum number of impurities in a cell to be 
<4) gives a total of 125 possible cell contents. The extra im- 
purities could be other types of atoms, but they could also be 



host vacancies V and interstitials /, although this would reduce 
the number of possible cell contents to 1 12 because introduc- 
ing a host vacancy and interstitial is the same as not intro- 
ducing either. Searching over the "impurities" H/N/0///V in 
silicon, with a maximum number of impurities of <4 and the 
same computational parameters as in the present study would 
cost about three times as much as the present study. For each 
set of cell contents we have in this work performed between 
about 20-900 structural relaxations, which appeared to be ad- 
equate in this case, although the required number is expected 
to increase rapidly with the size of the "hole" in the host crys- 
tal. In some cases it may be necessary to use larger "holes" to 
obtain all of the low-energy defects, and perhaps larger simu- 
lation cells as well. We have automated the search procedure 
and the calculations are performed in parallel. It appears to us 
that defect structure searches of the type described in this pa- 
per with up to, say, five types of impurity are perfectly feasible 
with modern computing facilities. 
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